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Loop quantum cosmology (LQC) provides promising resolutions to the trans-Planckian issue and 
initial singularity arising in the inflationary models of general relativity. In general, due to different 
quantization approaches, LQC involves two types of quantum corrections, the holonomy and inverse- 
volume, to both of the cosmological background evolution and perturbations. In this paper, using 
the third-order uniform asymptotic approximations, we derive explicitly the observational quantities 
of the slow-roll inflation in the framework of LQC with these quantum corrections. We calculate 
the power spectra, spectral indices, and running of the spectral indices for both scalar and tensor 
perturbations, whereby the tensor-to-scalar ratio is obtained. We expand all the observables at 
the time when the inflationary mode crosses the Hubble horizon. As the upper error bounds for 
the uniform asymptotic approximation at the third-order are < 0.15%, these results represent the 
most accurate results obtained so far in the literature. It is also shown that with the inverse- 
volume corrections, both scalar and tensor spectra exhibit a deviation from the usual shape at 
large scales. Then, using the Planck, BAO and SN data we obtain new constraints on quantum 
gravitational effects from LQC corrections, and find that such effects could be within the detection 
of the forthcoming experiments. 


I. INTRODUCTION 


The inflationary cosmology provides the simplest and 
most elegant mechanism to produce the primordial den¬ 
sity perturbations and primordial gravitational waves 
(PGWs) [U [2]. The former grows to produce the large- 
scale structure (LSS) seen today in the universe, and 
meanwhile creates the cosmic microwave background 
(CMB) temperature anisotropy, which has been exten¬ 
sively probed by WMAP [5], PLANCK [3], and other 
CMB experiments as well as galaxy surveys. PGWs, 
on the other hand, produce not only a temperature 
anisotropy, but also a distinguishable signature in CMB 
polarization-the B-mode, which was once thought to be 
already observed by BICEP2 [5], although subsequent 
analysis of multi-frequency data from BICEP2/Keck and 
Planck Collaborations showed that the signals could be 
due to galactic dust [B] (see also 0), and further con¬ 
firmation is needed. These observations in the measure¬ 
ment of the power spectra and spectral indices, together 
with the forthcoming ones, provide unique opportunities 
for us to gain deep insight into the physics of the very 
early Universe. More importantly, they may provide a 
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unique window to explore quantum gravitational effects, 
which otherwise cannot be studied in the near future by 
any man-made terrestrial experiments. 

Inflation is very sensitive to Planckian physics [5]. 
In particular, in most inflationary scenarios, the energy 
scale of the inflationary fluctuations, which relates to the 
present observations, was not far from the Planck scale at 
the beginning of inflation. As in such high energy regime 
the usual classical general relativity and effective field 
theory are known to be broken, and it is widely expected 
a quantum theory of gravity could provide a complete de¬ 
scription of the early Universe. However, such a theory of 
quantum gravity has not been established yet, and only 
a few candidates exist. One of the promising approaches 
is loop quantum gravity. On the basis of this theory, loop 
quantum cosmology (LQC) was proposed, which offers a 
natural framework to address the trans-Planckian issue 
and initial singularity, arising in the inflation scenarios. 
In fact, in LQC, because of the quantum gravitational ef¬ 
fects deep inside the Planck scale, the big bang singular¬ 
ity is replaced by a big bounce mm- This remarkable 
feature has motivated a lot of interest to consider the 
underlying quantum geometry effects in the standard in¬ 
flationary scenario and their detectability m- 

In LQC, roughly speaking, there are two kinds of 
quantum gravitational corrections to the cosmological 
background and cosmological perturbations: the holon- 
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omy [IMH] and the inverse-volume |19H24) . The main 
consequence of holonomy corrections on the cosmolog¬ 
ical background is to replace the big bang singularity 
of the Friedmann-Lemaitre-Robertson-Walker (FLRW) 
universe by a big bounce. The cosmological scalar m, 
vector [15], and tensor perturbations [IMl] have been 
calculated explicitly with holonomy corrections (see also 
m)- Due to these corrections, the modifications of 
the algebra of constraints generically leads to anomalies 
[mds]. Recently, it has been shown that these anoma¬ 
lies can be removed by adjusting the form of the quan¬ 
tum corrections to the Hamiltonian constraint. This is 
achieved by adding suitable counter terms that vanish in 
the classical limit dSHn]. With these anomaly-free cos¬ 
mological perturbations, the dispersion relation of the 
inflationary mode function (of scalar and tensor 

perturbations) is modified by the quantum corrections 
to the form [T8| , 

( 1 . 1 ) 

where pc is the energy density at which the big bounce 
happens. In the classical limit p ^ pc, the above equa¬ 
tion reduces to the standard one of GR. From this equa¬ 
tion, one can see that when the energy density p ap¬ 
proaches Pc/2, the dynamics of cosmological perturba¬ 
tions is significantly modified by the holonomy correc¬ 
tions. With this modified dispersion relation, the power 
spectra for both scalar and gravitational wave pertur¬ 
bations were calculated up to the first-order approxima¬ 
tions of the slow-roll parameters, and the observability of 
these corresponding quantum gravitational effects were 
discussed in some detail [nmsi. 

The inverse-volume corrections are due to terms in 
the Hamiltonian constraint which cannot be quantized 
directly but only after being re-expressed as a Poisson 
bracket. It was demonstrated that the algebra of cos¬ 
mological scalar [191 EO], vector [21], and tensor per¬ 
turbations |22j with quantum corrections can be closed. 
Consequently, the cosmological perturbed equations for 
scalar and tensor perturbations are modified [531 HI]. For 
the scalar perturbations, the dispersion relation takes the 
form [53] 



( 1 . 2 ) 

where the constants ao, do, and a encode the specific fea¬ 
tures of the model, and Spi{r]) is time-dependent, which 
usually behaves like dpi ~ a~'^. For tensor perturbations, 
ojlip) is given by [53], 

Wfc(?7) = (l + 2aodpi(77))^^- (1-3) 

The power spectra of both scalar and tensor perturba¬ 
tions due to the inverse-volume corrections, again up to 


the first-order approximations of the slow-roll parame¬ 
ters, were studied in [53], in which some constraints on 
some parameters of the model were obtained from obser¬ 
vational data [53]. The non-Gaussianities with inverse- 
volume corrections was also discussed in [53] . 

Although a lot of effort has already been devoted to 
the inflationary models of LQC with both holonomy and 
inverse-volume quantum corrections, very accurate cal¬ 
culations of inflationary observables in LQC are still ab¬ 
sent, and with the arrival of the era of precision cosmol¬ 
ogy, such calculations are highly demanded. Recently, we 
have developed a powerful method, the uniform asymp¬ 
totic approximation method [57IlH()j . to make precise ob¬ 
servational predictions from inflation models, after quan¬ 
tum gravitational effects are taken into account. We 
note here that such method was first applied to inflation¬ 
ary cosmology in [3T], and then we have developed it to 
more general mult-turning points cases and more precise 
higher order approximations |27H30j . The main purpose 
of the present paper is to use this powerful method to de¬ 
rive the inflationary observables in LQC with holonomy 
and inverse-volume quantum corrections with high accu¬ 
racy. More specifically, we consider the slow-roll inflation 
with the quantum gravitational corrections, but ignore 
the pre-inflation dynamics. By using the general expres¬ 
sions of power spectra, spectral indices, and running of 
spectral indices obtained in [53], we calculate explicitly 
these quantities, up to the third-order approximations 
in terms of the uniform asymptotic approximation pa¬ 
rameter, for which the upper error bounds are 0.15%. 
All the inflationary observables are expressed in terms of 
slow-roll parameters and parameters representing quan¬ 
tum corrections explicitly. In the present paper, we also 
provide the expansion of all these inflationary observables 
at the time when the inflationary mode crosses the Hub¬ 
ble horizon, and calculate the tensor-to-scalar ratio. It is 
interesting to note that the holonomy corrections do not 
contribute to the tensor-to-scalar ratio up to the third- 
order approximation. More interestingly, it is shown that 
with the inverse-volume corrections, both scalar and ten¬ 
sor spectra exhibit a deviation from the standard one 
at large scales, which could provide a smoking gun for 
further observations. 

The paper is organized as follows. In Sec. H, we 
present all the background evolution and perturbation 
equations with the holonomy corrections, and calculate 
explicitly the power spectra, spectral indices, and run¬ 
ning of spectral indices. Then, in Sec. HI we turn to 
consider both background and perturbations with the 
inverse-volume corrections, and calculate all the infla¬ 
tionary observables. Our main conclusions are summa¬ 
rized in Sec. IV. Three appendices are also included. In 
Appendix A, we give a brief introduction to the uniform 
asymptotic approximation method with high-order cor¬ 
rections, while in Appendices B and C, we present some 
quantities discussed in the content of the paper. 

Part of the results to be presented in this paper was 
reported recently in [32] . In this paper, we shall pro- 
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vide detailed derivations of these results, and meanwhile 
report our studies of other aspects of LQC inflationary 
cosmology. 


II. INFLATIONARY OBSERVABLES WITH 
HOLONOMY CORRECTIONS 

In this section, let us consider the inflationary cosmol¬ 
ogy with the holonomy corrections. 


A. Background equations and equations of motion 
for scalar and tensor perturbations 


To begin with, let us first consider a flat FLRW back¬ 
ground 


ds^ = a'^{ri){—dri^ + dx^dx^), ( 2 . 1 ) 


where 0 ( 77 ) is the expansion factor and 77 the conformal 
time. With the holonomy corrections, the Friedmann 
equation is modified to the form m 





( 2 . 2 ) 


where H = a/a is the Hubble parameter with a dot 
representing derivative with respect to the cosmic time 
t {dt = adrj), p is the energy density of the matter 
content, and pc is a characteristic energy scale of the 
holonomy corrections and usually is of the order of the 
Planck energy density: Pc TTZpj with the Planck mass 
TOpi = 1.22 X 10^®GeV. A general prediction associated 
with the above equation is that the holonomy corrections 
lead to a resolution of the big bang singularity, in which 
it is replaced by a non-singular big bounce occurring at 
p = Pc- For a scalar field ip with potential Vipp), the 
Klein-Gordon equation reads 


ii + 3ir^ + tW)=„, 

dip 

The energy density of the inflaton field is 

■2 


P = 




■Viif). 


(2.3) 


(2.4) 


The above set of equations determines uniquely the 
evolution of the FLRW background. As shown in [HH] . 
“the standard slow-roll inflation” can be triggered by the 
preceding phase of quantum bounce. In this paper, for 
the sake of simplification, we shall focus on the slow-roll 
inflation with the holonomy corrections, and ignore the 
pre-inflation dynamics. In this case, we are in a process 
such that the energy density of the cosmological fluid is 
supposed to be dominated by the potential of inflaton ip, 
i.e., ip^ <C V{ip). With this condition, it is convenient to 
define a hierarchy of Hubble flow parameters. 


eo = 


H 


^n+1 


din e„ 
din a 


(2.5) 


On the other hand, as the quantum holonomy correction 
is also small, it is convenient to introduce the parameter 
by, 

6h = — <^ 1. (2.6) 

Pc 

Then we have 0=1 — 2p/pc = 1 — 2dpr. 

Let us now consider cosmological perturbations. With 
the holonomy corrections, the equations for cosmological 
scalar perturbations can be cast in a modified Mukhanov 
equation for the gauge-invariant mode function 

[Min], 


d^Pk\v) 

dr]^ 



zM) 

Zsiv)/ 


P-kHv) = 0, 


(2.7) 


where 07 ^( 77 ) = rt{r])k^, and Zs = aip/H. Similarly, the 
equation for the tensor perturbation can be cast in the 
form [17], 


drj'z 




where Zt = a/VTi. 


B. Power spectra and spectral indices in the 
uniform asymptotic approximation 


To apply for the uniform asymptotic approximation, 
first we write the equations of motion for both scalar (Eq. 
(2.7)) and tensor (Eq. (2.8)) perturbations in the form 
(A.l) by introducing a new variable y = —kp. Then, we 


find that g{y) and q{y) must be chosen as [571 HE], 




z"{y) 


1 


zirf) J 47/2 


12 ^ {p) 


qiy) = - 


yz 

1 

4y2’ 


-clip), 


(2.9) 

( 2 . 10 ) 


where Cs{p) = \/^ip) and v^{p) = p"^z”(p)/z{p) + 1/4. 
Eor the scalar pertur bati ons, u sing the slow-roll parame¬ 
ters defined in Eqs. (2.5) and (2.6), we find that 


f(d) 




Zs{r]) 


€263 

2 


ei£2 

2 


— GeiSn + 6eid// — 4eie2(5/7 — 18eid|^^. 

( 2 . 11 ) 

In the above, we used the subscript (or superscript) “s” 
to denote quantities associated with the scalar perturba¬ 
tions. Similarly, for the tensor perturbations, we find 


zj'iv) 

ztiv) 


~ a^H^(2 - ei - eSnei - ISSj^ei - 
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+ 38(5 


2 2 
H^l 


— 25 //eie 2 — 6 ( 5 |^eie 2 ^ , 
( 2 . 12 ) 


lim /a = -vq 
y-s-O 


where the subscript (or superscript) “t” denotes quanti¬ 
ties associated with the tensor perturbations. 

With the functions g{y) and q[y) given above, we are in 
the position to calculate the power spect ra an d spectral 
i ndice s from the general formulas Eq. (A. 26 ) and Eq. 
( A. 27 |). As we discussed in [ 55 ], in order to do so, first we 


- 12 In 2 
24 


-^0 1 1 - ^ - ln2-h 


Cl _ 
Co 

ln^2 


f'o 

C2 

Co 


need to expand v{rj) and 0^(77) around the turning point 
yo of g{y) (i.e., g{ yo) = 0 with = - fc??o), then perform 
the integral of \/g{y) in Eq. (A. 26 ), and calculate the 
error control function J^(+oo). In the slow-roll inflation, 
it is convenient to consider the following expansions, 


f C( 3 ) 7 r 2 ln 2 ln ^3 1 , 3 y ^ _ 

+ a)"- 

(2.18) 

In the above, C(^) denotes the Riemann zeta function. 

Now, we turn to consider the error control function 
which in general can be written as 


2 2 - , - 1 2/ ,1- 1 2 y 

v{ri) ~ 12q-I- 123 In —-I--122 In —, 

yo 2 yo 

2 X - , - 1 y ,1- 1 2 y 

Cs(y) Ci Co-I-Cl In —-I--C2 In —, 

yo 2 yo 


( 2 . 13 ) 

( 2 . 14 ) 


mo = 


yo 


'yo 


5^/2 


16 ^' 


,3 + ^ f \^ 9 dy . (2.19) 


with 


In the limit y —)■ 0 , after some lengthy calculations, the 
above expression can be cast in the form 


Vl = 


dv{g) 


dln(—77) 



d^i^(g) 


c;2^(+Oo) 


7 

^0 

dln^(—77) 

7 

^0 

A 

6720 V Co y 


i>i (23 + 121 n 2 ) 


721^2 


dcs( 77 ) 

Cl = 


dln(—77) 


_ _ d?Cs{g) 
C2 = 


din (—77) 


( 2 . 20 ) 


(2.15) 


In the above, 770 is the conformal time for mode k at 
the turning point yg. In the slow-roll inflation with the 
holonomy corrections, in general we have 12(77) ~ ^+ 0 {e), 
Di ~ 0 {e‘^), V2 — and co ~ 1 -|- 0 {Sh), ci ~ 

0 {€Sh), C2 — 0 {e^ 6 H)- The slow-roll expansions of all 
these quantities are presented in Appendix C. 

With the above expansions, we notice that \/g{y) = 
yj\^g{y) can be expanded as 


As we only need to calculate the power spectra up to 
the second-order in slow-roll parameters, in the above 
we have ignored cf, C2,72^, 722, terms. Once we get the 
integral of yjg(y) (as presented in Eq. ( 2 . 17 )) an d err or 
control function Jf'(-l-oo) (as presented in Eq. ( 2 . 20 )), 


from Eq. (A. 26 ) we can calculate the power spectra. 


To calculate the corresponding spectral indices, let us 
first consider the /c-dependence of 729(770), 721(770) through 
do = Voik)- From the relation —kijQ = 720(770)700(770), we 
observe that 


\/5(y) - 




+ 


729721 


yy/vl-cly^ 


cociy^ In — 

yo 


dln(- 77 o) 
din fc 


~ -1-h 


Cl 

721 

r'l _ 


Co 

720 

Vco 

720 / 


■ ( 2 . 21 ) 


720722 


CoC2y 


Then, the spectral index is given by 


2yyjvl - cly'^ 2 1/122 _ ^2y2 


(2.16) 

Therefore, the integral f y/gdy is divided into three parts, 
rvo 

/ V 9 iy)dy = h+h +I3, ( 2 . 17 ) 

^ V 


n-1 ~ (3-27 >o)+^^+ f77^-21n2)7?i 


1 


2720 


Co 6 co 72 o 


672, 

2Po In 2 
Co 


C2 


23 + 121n2 , TT^ ,2 


where 




Similarly, after some tedious calculations, we find that 
the running of the spectral index a = dn/dlnk can be 
written in the form 


lim Ji 
y-s-o 


lim I2 
y-s-o 


(1 — ln2)ci72o / TT^ 


ln^2 


2 yo 


^ 1 , 


2 vocl 4721C1 

Q!(fc) ~ -z- 


'.'0 


Co 

2 In 2720 2720 


1 2.2 

3 P 3 ^--- 


Co 


'^0 ^0 
1 

Co ' 6coPo 


C3 + 2721 
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1/1 1 \ - 2 C 2 J/o 


Co 


ln^2 


TT^ ln2 23 \ 


12 




Vs- 


(2.23) 


C. Scalar perturbations 

With the slow-roll expansions of Pq) vi, i> 2 , and 
Cl, C 2 , C 3 presented in the above and after some tedious 
calculations, we obtain the scalar spectrum, 

(^) — "I ^ ~ ^ + 1 ) Cl ~ + 2 

— 2 ( 2 Z)p + 3 ) — Dp 6 H ^2 

+ ^2l)p + 2Z)p + ^— 5 + Ai^ ef 
+ f-Dp — zip + — 8 + Ai + 2 A 2 ) €162 


2^^ 


TT^ 3 Ai \ _2 


7r2 - iZ 2 \ 

-+Ai-— 62^3 


(2.24) 


where Ag = 


— ISIH^ 


72e37r2ei > 


- ln2, Ai = 1®3606 


181 


32761 


and A 2 = 


9269 


Note that a letter with an over 


2 ’ “ ^ “ 589698 ■ 

bar denotes a quantity evaluated at the turning point j/o- 
Then, the scalar spectral index is 

— 1 — 2ei — e2 + — 2ef — Dn£2£3 

— (2Dj^ + 3 ) 6162 + 12 ^|fei + AD^5 — 26 ^ 

—2 (Z>n — 1 ) 5//eie2 


-{6Dn- 2A„1 - TT^ + ^ 


e?e 2 


- I 3iZ„ 

Anl 

~Y~ 

- ( 


^"-^-A„i- 2 A„ 


25 

y 


£162 


TT 

T 


eies 


y"-^-A„i- 2 A„ 


22 

n 


An2-^ + C2e3e4), 


^ eie2e3 


(2.25) 


where Z)„ = ^ - ln2, A„i = and A „2 = 

The running of the scalar spectral index reads 

as ~ — 2 eie 2 — £263 + 4ey_f7 — 2 ei£2^/7 — 66^2 

— (3 + 2Dn) €i £2 — 2 (z)n + 2) £l€2£3 ~ 

—Z)ne2e3e4 — SiD^elSn — i2ele2 


+ (27r2 - 14£»„ + 4 A„i - 39) efel 

'lir^ - - - - 25 

- 3Dn — + Anl + 2An2-^ 

65 \ 

+ ( 77^ — 8Z)n + 2Ani- — j £^2^3 

+ ( — - 105„ - 3Dl 

71 \ 

+3Ani + 6An2-j ClC^CS 

Ttt^ 22' 

-5Zin ~ + Anl + 2An2-) eie2C3 


t2 .D? 


+ ( ^ + Anl - 5j £^£3 + 1 ^—-^ + An 2 ) £263 

+ ( -5Zin — + Anl + 2An2-) eie2C3C4 

_L I ^ _L Am 8 \ 2_ _ 

+ II + — 


2 3D^ 


+ -^ + 3An2j £2e3£4 

+ -^ + An 2 ^ £26364 + 6 (2Dn + 3 ) £^2^77 

+ f ^-^ + An2 j e2e3£4e5 ~ 2 (Z)n + 3 ) €i£2577 

— 2 (zin + 2) £l£2e3^7/ + — 6£i£2^77. 

(2.26) 

D. Tensor perturbations 


Similar to the scalar perturbations, within the slow- 
roll approximations, we find that the power spectrum for 
the tensor perturbations reads 

At(^) — / 4 t|l -|-^77 — 2 (iZp-|-l)ei 

3^2 

H—^ — 2{2Dp -j- 3)SHii 


Ai + — -5 + 2Dp + 2D^p] £f 


2 A 2 - 2 - 2Dp - 




(2.27) 


where At = Also the tensor spectral index and 

its running are given by 

Tit — —2£i -|- 4577ei — 2ei — 2[D^ -I- l)eie2 + 125|i£i 
-l-dZin^Ttei — 2 ei — 2 (L>n — 1)^776162 
7 - - 50 \ 

TT^ — 6Z)n + 2Anl-) 6^2 


— — 2Dn — Z)n + 2An2 — 2 £i £2 
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12 


— 2 _Dn — + 2 An 2 — 2 ] £ 16263 , 


(2.28) 


and 


at — —26162 + 4 ( 5 //ei — 2d^reie2 — 66162 — 2(i)„ + 1)6162 
—2(Z)n + 1)616263 ~ SiDnSn^i + 2(7 + 6lin)feeie2 

—2{Dn + 2 )^ 1/6162 — 2(i)n + 2)d//eie2e3 — 12 eie 2 


+ (27r2 - 36 - UDn + 4Ani) 
^2 


+ 


TT 


- 2 L»n - £»„ 


12 

+ ( TT^ — szin 4“ 2 Anl — 


2An2 ~ 2 I £ 16 : 
56 


eie2£3 


TT 


TT 


— 6i4„ — 3I1„ + 6An2 — 6 ) £i£2£3 


12 

TT^ 

I 2 


— 2i4„ — + 2A„2 — 2 


eie2e3 


— 2Dn — + 2An2 — 2 I £i62£3£4- (2.29) 


obtains 

Us — 1 — 26*1 — 6*2 — 26*]^ — (3 + 2Z1*)£*i£*2 

— 71n£*26*3 + 4(5 *h 6*1 + 12d*H£*l 

+ - 6 L>* + TT^ + 2 A* 1 ^ e*ie *2 - 2 £*i 

+ (^3^ ^^n2- ^ -311* - H*^^ 6 * 16*2 

+ ■*" ^^n2-^-411* - H*^^ 6*16*26*3 

+ 4—y “ 2^ £* 2^*3 4- (l2Zin - 811*) (5 *h£*i 

/tT^ £)*2\ 

4- ( y 4- A *2 -^ j (e* 2£*3 4- £* 26 * 36 * 4 ) 

+ { AD * - 65 n + 2 ) ( 5 *h 6 *i£* 2 . ( 2 . 31 ) 

The running of the scalar spectral index reads 

Cts — —26*16*2 — 6 * 26*3 — 2(Z1* + 2)£*i£*2£*3 

— 6 £*i 6*2 — I 1 *£* 2 £* 36*4 — (3 4 " 2 Il*)£*i £*2 

— Z1*£*26*3 4- 4i5*h£*i — 2i5*h6*i£*2 — 126*i£*2 


E. Expansions at horizon crossing 


In the last two subsections, we have obtained the ex¬ 
pressions of the power spectra, spectral indices, and run¬ 
ning of spectral indices for both scalar and tensor per¬ 
turbations. It should be noted that all these expressions 
were evaluated at the turning point yo- However, in the 
usual treatments, all expressions were expanded at the 
horizon crossing a{r]^,)H{ri^,) = Thus, it is use¬ 

ful to rewrite all the expressions given in the last section 
at the time when the scalar (or tensor) mode crosses the 
horizon. After some tedious calculations, for the scalar 
perturbations, we find the scalar spectrum can be ex¬ 
pressed as 


A 2 (fc) cs a: 1 - 2 (1 + H;) 6*1 - zz;e *2 + d*H 


2D;^ + 274; + - - 5 + AC ) e: 

1 2 \ 2 

2 < + y-' + xh- 


4 " 2^*^ ~ 44 ;d*H £*2 — (471; 4 - 6) i 5 *h 6 *i 
4 - {— D * + —- 7 4 - a; -I- 2 A;^ £*i£*2 


1 


y “ 2^^? ^2 ) 6*26*3 


(2.30) 


where the subscript denotes evaluation at the horizon 

crossinp- A * = n* _ _ 6 L _ In 3 A* — "^^5296 _ 

crossing, _ Y 2 e 37 r 2 e*i I “ 181 ^1 “ 98283 

y, and A 2 = gggggg . Foi thc scalar spectral index, one 


4 - ^ 2 ^ -I- 4 A;;^ ^ 4474;^ 6*16*2 

+ + AOi + 2 a ;2 

-|-374;-74f^ 6*1632 
4 - -I- 2 A;^ —-— 874 ;^ 6*16*26*3 

4 - + 3 A;i 4 - 6A;2 - 23 


— 1074; — 374;^ 6*16*26*3 — 6(5*ji6*i£*2 


77r^ 

23 


AC, + 2 A ;2 - 574 ^ - D * 


A* 


A* — 

^n 2 


6*16*26*3 

TT^ 14 

y ~ y 


74 


*2 


TT 

y 


£2 f2 

6 * 26*3 


6*26*3 


7tt^ 

u 


A*i + 2A 


n2 


-^- 574 ;- 74 ; 


6 * 16 * 26 * 36*4 


A: 


nl 


7 


2 y “ 31 ^* 26 * 36*4 4 - i6d*H£*i 


3 --,*2 


3a;2 - 2 ^^^ ~ ^ 6 * 26 * 36*4 
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+ 2 ^ (e*2e*3e*4 + e*2e*3e*4e*5) 

—2 {D* + 3) i5*He*ie*2 ~ ^ {D* + 2) (y*He*ie*2e*3 
-8L»*5*He*i + 6 (27?* + 3) 5,He*ie*2. (2.32) 

Similar to the scalar perturbations, now let us turn to 
consider the tensor perturbations, which yield 

A^{k) ~ 11 — 2 (l + i?*) e*i + (5 *h + 2^*h 

+ (^2 d;’^ + 2d; + ^-5 + a{^ 4 

+ “ 2i?* + 2 + 2A2^ e*ie*2 

-2(27?; + 3)^*He*i|. (2.33) 

For the tensor spectral index, we find 

Tit — ~2e*i — 2(1 + i?*)e*ie*2 ~ 2eJ]^ + 4(5*He*i 
~2e*i + — 6D* + 2A*;^ —e*ie*2 

+ (^^ - 2 + 2A*2 - 27?: - e*i4 

+ ^Y2 ~ ^ 2A:2 — 27?: ~ e*ie*2e*3 

+ (124 - 87?:) 

+ (47?: ~ 67?n + 2) <5*HC*ie*2 + 12^4£*i- 

(2.34) 

Then, the running of the tensor spectral index reads 

Oit — ~2e*ie*2 — 2(1 + 7?:)e*ie*2e*3 ~ 6eJj^e*2 
~12e*4e*2 ~ 2(1 + 7?:)e*ieJ2 + 4(5*He*i 
~2(5*He*ie*2 + 16<54£*i ~ 87?:(5*He*i 

+ ^27r^ — 147?: + 4 A:j e*ie*2 

+ (^^ + 2A:2 - 2 - 27?: - 4") 

+ “ 87?: + 2A:;^ —e*ie*2e*3 

+ — 6 + 6A:2 — 67?: — 37?:^^ e*ie*2£*3 

+ ^Y2 “ ^ “ 27?: “ 4^^ e*ie*24 

+ ^Y2 ~ ^ 4” 2A:2 — 27?: ~ 4^^ e*ie*2e*3e*4 
T 2 (67?+ 7) i5*hc*i£*2 

—2 (7?: + 2) (5*He*ie*2 ~ 2 (7?: + 2) (y*He*ie*2e*3. 

(2.35) 


Finally with both scalar and tensor spectra given 
above, we can evaluate the tensor-to-scalar ratio at the 
horizon crossing time 77 *, and find that 


r 


16e*i 


1 + 4^*2 + 




£162 


+ 



e 


2 

*2 



(2.36) 


It is remarkable to note that the holonomy correction 
parameter <5*// doesn’t contribute to the tensor-to-scalar 
ratio r, up to the third-order uniform asymptotic approx¬ 
imation. 

In addition, to the first-order of the slow-roll param¬ 
eters, it can be shown that our results given above are 
consistent with those presented in . 


III. INFLATIONARY OBSERVABLES WITH 
INVERSE-VOLUME CORRECTIONS 

Now let us turn to consider another type of quantum 
gravitational correction, the inverse-volume, in LQC. 


A. Background evolution and equations for 
perturbations 

In the presence of the inverse-volume corrections, the 
effective Friedmann and Klein-Gordon equations read 

m 


= 


SttG 


4 

27? 


V{ip) 


rr,n ^dlni}\ 


d(f 


(3.1) 

(3.2) 


with p = and 


a ~ 1-I-ao<5pi, i? ~ 1-I-79o<5pi, (3.3) 


where (5pi characterizes the inverse-volume corrections in 
loop quantum cosmology and 

<5p,^(^) . (3.4) 

Note that here we only consider the inverse-volume cor¬ 
rection i5pi at the first-order C?((5pi). Thus to be consis¬ 
tent, through the whole paper, we shall expand all the 
quantities at the first-order of ^pj. In the above, oq, 
i?o, and opi are constants and depend on the specific 
models and parametrization of the loop quantization. 





Specifically for the parameter a, different parametriza- 
tion schemes shall provide different ranges of a [23l |34] . 
Moreover, oq and -do are related by the consistency con¬ 
dition 


i9o(cr — 3)(a -I- 6 ) = 3ao(cr — 6 ), 


(3.5) 


while a takes values in the range 0 < ct < 6 . 

The evolution of the background, which can be deter¬ 
mined by the above set of equations, is usually differ¬ 
ent from the evolution given in the standard slow-roll 
inflation, because of the purely geometric effects of the 
inverse-volume corrections. However, as indicated in |24| . 
in that regime, the constraint algebra has not been shown 
to be closed. One way to consider the slow-roll inflation 
with inverse-volume corrections is in the large-volume 
regime, where the quantum corrections are small and the 
constraint algebra is closed. In this paper, we will focus 
on the latter case. Similar to the case with the holon- 
omy corrections, we still adopt the slow-roll parameters 
defined in Eq. (2.5). 


The inverse-volume corrections can be also introduced 
into equations governing the evolution of cosmologi¬ 
cal perturbations. In particular, it was found that, 
when inverse-volume corrections are present, the gauge- 
invariant comoving curvature perturbation TZ is con¬ 
served at large scales. Such a feature of TZ strongly sug¬ 
gests that one can write a simple Mukhanov equation in 
the variable (ry) = ZgTZ, which is [521 133 


dif 


+ s^{r])k^ - 


Zs{v) 


(3.6) 


where 


Zsiv) 


aip 

~H 



otQ — 2i9o 
2 



(3.7) 


depends on the evolution of the background and 


s^(r/) = 1 + xdpi, (3.8) 


with 



For the cosmological tensor perturbations hk, when the 
inverse-volume corrections are present, the corresponding 
Mukhanov equation for the variable = Zth^. is 

written as [55] 


drf 


= 0, (3.10) 


with a{ri) being given in Eq. (3.3) and 


Zt 




(3.11) 


B. Power spectra and spectral indices in the 
uniform asymptotic approximation 


Similar to the last section, to apply the uniform asymp¬ 
totic approximation, we first write the equations of mo¬ 
tion for both scalar (Eq. (3.6)) and tensor perturbations 
(Eq. ( |3.10| )) into the standard form Eq. ( |A.1[ ). Then, 
for the scalar perturbations, the functions g{y) and q{y) 
must chosen as mm, 


A"5(y) = 

-lA 

s‘^{g)k‘^ — 

Zsiv)) 

1 

4y2’ 

(3.12) 

qiy) = 

1 

4y2 

5 



(3.13) 

for tensor perturbations one 

chooses 



A^5(y) = 

-h( 

a^{r])k^ — 

znv)\ , 

ztiv) ) 

1 

4y2’ 

(3.14) 

qiy) = 

1 

4y2- 




(3.15) 


In the slow-roll approximation, z'^/zg and z'ljzt can be 
casted in the form, 


z!!M 

Zsiv) 




3e2 


2-ei 


£162 , e| , £263 
2 4 2 


(3.16) 


and 

where 


4'iv) 




_ cr^(cr — 3 )q;o 


2 - £1 +/(‘H£,)<5pi1 , (3.17) 


-b 


4£i 

a{(j — 3)ao £2 


a{(j - 3)(cr -I- 6 )^?o 
12 


(T^ao 


£1 




cr(cr — 3) 


CkQ. 


(3.18) 


Because 5pi ^ it is convenient to write the function 
g{y) in a simplified form 


^ 9\y) = —5 - 1 - X<JP 1 H- 3 —Opi, 


y. 


y" 


(3.19) 


where in a slow-roll background, v{rj) and m{ri) are slow- 
rolling variables depending on the types of the pertur¬ 
bations. In particular, for the scalar perturbation, y is 
given by Eq. (3.9), while for the tensor perturbation we 
shall replace y by 2 ao. 

With the functions g{y) and q{y) given in the above, 
we are in the position to calculate the power spectra and 
spectral indices from the general formulas Eq. (A.26) and 
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Eq. (A.27|. However, unlike the case for the holonomy 


corrections, in which i/(ry) and ^( 77 ) are slow-rolling quan¬ 
tities, ^pi in Eq. (|3.19| cannot be treated as a slow-rolling 


qua ntity duri ng infla tion. Conse quently, the expansion 
of \/ giy) = in Eq.(2.16 ) can not be directly ap¬ 

plied to the function g{y) of Eq.(3.19|. 


In order to apply for the formulas Eq.(A.26) and 


Eq.(A.27) to calculate inflationary observables with the 


inverse-volume corrections, let us first write 5pi as 


5p. = 


H 


V 7. y 

\ k J \—ar]HJ 


y'^ = epiK(? 7 )y‘^, (3.20) 


with epi = {^Y ^ ^ ^od assume 

that tr is an integer within the range 0 < cr < 6. With 
these conditions we can write the function g(y) in the 
following form 


giy) = 


yo-y 

y2 


(ho + hiy + 


+ h^y'^ + h^+iy'^'^^) . 

(3.21) 


In the above yo = yoiv) is assumed to b e slow-rolling. 
Comparing the above form with Eq. (3.19) we find 


h(j+i 

hfj 

h(7— 1 


= Xepi«, 

= xyoepiK, 

= ixyo - myo^iyoepiK, 


hi = (x2/o - "1% )yo (cr>j>2). 


h2 = ixyo - myQ^)yQ~‘^epiK, 
hi = 1 + ixyo - myQ^)yQ~^epiK, 
ho = yo +ixyo-myQ^)yQepiK. 


(3.22) 


Note that in the above we have assumed that m = 771 ( 77 ) 
and K = nit]) are slow-rolling quantities. Then, expand¬ 
ing Ygiy) in terms of epi to the order C>(epi), we find 




epiK 


o xivY^ + yoy 

wyo + V L 


+ --- + y^n 


-"^ivY^+Voy 

ho 


,,<^+1 


(3.23) 


and 



(3.24) 


yY^ + yoy 


= yo""^ 


+ yly"~^ 


, u a I s 

a + 0 -he — 


yo 


yo 


+ e 





3 


(3.25) 


where the relations between the values of cr and 


{a, 6 , c,(I, e,/} is 


a = 6 


10 

II 

b 


b 

II 


b 

II 


b 

II 

to 


(7=1 



{ 1 , 1 , 1 , 1 , 1 , 1 }, 
11 , 1 , 1 , 1 , 1 , 0 }, 
11 , 1 , 1 , 1 , 0 , 0 }, 
11 , 1 , 1 , 0 , 0 , 0 }, 
11 , 1 , 0 , 0 , 0 , 0 }, 
11,0,0,0,0,0}. (3.26) 


Note that when a is an integer, one can always set a = 1. 

Now let us turn to the integral of \/giy)- In order 
to carry out the integration, we need also to specify the 
form of all the slow-rolling quantities yoiv), n^iv), and 
K{r]). To do so, similar to the case with the holonomy 
corrections, it is convenient to expand all these quantities 
around the turning point yo of giy) (i.e., giyo) = 0 with 
yoivo) = —krjo) in the slow-roll inflation, i.e.. 


yoig) 

- 2/0 + yiln—, 
yo 

(3.27) 

777 ( 77 ) 

~ 7770 + In —, 

yo 

(3.28) 

7C(77) 

~ Ko -f 7 ti In —. 

(3.29) 


yo 


In general the quantities yo, yi, mo, mi, ko, ni are of 
order of 


yo - 

- 2 + 0 (e 



yi - 

- oiY) + 

0 (epi), 


7770 ^ 

- Oij), 

7771 ^ 0(1), 


Ko - 

- HOil), 

Ki ~ HOiti). 

(3.30) 


As we have assumed 0 < ct < 6 , the two sequences in the 
above expressions are finite and can be expressed as 


+ yo y + • • ■ 

+ y"+' 


= yo^^ 

i-f 

y , 

-ha 

(- 


yo 

\yo 


- 1 - c 

(s) 

+ d 




5 


As shown in [24) . when cr is in the range of 1 < cr < 6 , 
the parameter of the inverse-volume corrections i7*epi 
is constrained by the observations that i7*epi should be 
< 10“^. As we shall show below, in our calculations, the 
constraint is tighter because of enhancement of the 
inverse-volume corrections. In this situation, as the slow- 
roll parameter is usually at order of 10 “^, the correction 
term H'^epi is expected to be < ef, which usually is at 
the order of 10“^. Thus, in this paper we consider the 
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inverse-volume corrections H'^epi as the second-order in 
the slow-roll expansion. 

Then, with these expansions we can perform the inte¬ 
gral of Eq. (A.26 1 , and calculate the error control func¬ 
tion ^{+oo). The results of the integral of \/g{y) and 
the error control function J^f(+oo) are all presented in 
Appendix B. The slow-roll expansions of toq, toi, and 
Ko, Ki etc are presented in Appendix C. 


C. Inflationary spectra for both scalar and tensor 
perturbations 


+ ^^1 + -^— 5-1- 2Dp + 2Ilp^ €i 
+ ^ 2 A 2 ~ ^ ^2 ~ ~ ^ 1^2 

+ epi0ij) (3.34) 

Tit — ~2ei — 2eJ — 2 -|- l) eie2 + epi 

(3.35) 


Let us first consider the scalar spectrum. With the 
slow-roll expansions of Dq, ui, V 2 , mo, mi, m 2 , etc, 
present ed in Appendix C and the expression of yo given 
in Eq. (B.7), we find that up to the second-order in the 
slow-roll parameters, the scalar spectrum can be cast in 
the form 

Alik) ~ A J 1 - 2(1 + .Dp)ei - .Dpe2 


2D„ 


2Zl^+--5 


+ I 


7tt^ 
2 3 


-8-h Ai -H2Ai 


-L )2 
2 P 

TT ^ -r 1 - 2 

^ + 2 - 2 ^p 


TT* 6 1 . 

¥-2 + 3^' 


£263 


epi 


3H 

~Y 


ei 


eie 2 




fell £2 


£1 


(3.31) 


Then, the corresponding scalar spectral index and run¬ 
ning read 

u-s — 1 ~ 2ei — £2 ~ 2 e^ — {2Dn + 3) eie 2 ~ DnC2£3 


(3.32) 



as — —2eie2 — £2£3 



(3.33) 


where Q^i, 


Q, 


0 ’ 


and 


are given in Table II. 

Similar to the scalar perturbations, up to the second- 
order of the slow-roll parameters, the tensor spectrum, 
spectral index, and the running, are given, respectively, 

by 


Alik) ~ AJ 1-2(1 + .Dp) ei 


and 


at 


—2eie2 -I- £pi 



(3.36) 


where Qq\ and £q^ are given in Table II. 


D. Evaluating at horizon crossing 


So far, we have obtained all the expressions of the 
power spectra, spectral indices, and running of spectral 
indices for both scalar and tensor perturbations with the 
inverse-volume corrections. However, to compare with 
observations, we need to express them in terms of the 
slow-roll parameters which are evaluated at the time 77 * 
when the scalar or tensor modes cross the horizon, i.e., 
a(? 7 *)iL(r 7 *) = s(? 7 *)A: for the scalar perturbations and 
a(? 7 *)iL(r 7 *) = 0 ( 77 *)^: for the tensor perturbations. Be¬ 
cause in general the values of s(? 7 *) and 0 ( 77 *) are dif¬ 
ferent, for the scalar and tensor modes with the same 
wavenumber k, they cross the horizon at different times. 
When s(? 7 *) > 0 ( 77 *), the scalar mode leaves the hori¬ 
zon later than the tensor mode, and for 5 ( 77 *) < a(? 7 *), 
the scalar mode leaves the horizon earlier than the ten¬ 
sor. In this case, as we have pointed out in m , caution 
must be taken for the evaluation time for all the infla¬ 
tionary observables. As we have two different horizon 
crossing times, it is reasonable to re-write all expressions 
in terms of quantities evaluated at the later time, i.e., 
we should evaluate all expressions at scalar-mode hori¬ 
zon crossing a(? 7 *)iL( 77 *) = 5 ( 77 *)^ for 5 ( 77 *) > a(? 7 *) and 
at tensor-mode horizon crossing airj^,)Hirj^) = air]^)k 
for 5 ( 77 *) < 0 ( 77 *). However, detailed analysis shows that 
such a difference only contributes to the high-order terms 
in terms of the slow-roll parameters, which is beyond the 
approximation we consider here. Thus, in this paper, 
we will not distinguish these two different cases and only 
consider the expansions at the time when the scalar mode 
crosses the Hubble horizon. 

Then, we shall re-write all the expressions in terms of 
quantities evaluated at the time when the scalar mode 
leaves the Hubble horizon a( 77 *)iL(? 7 *) = s( 77 *)fc. Skip¬ 
ping all the tedious calculations, we find that the scalar 
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spectrum can be written in the form 


^*(S) ^ 


(3.42) 


Alik) A:\l-2{1 + D;)e,,-D;e,^ 


For the tensor spectrum, we get 


2 Dl^ + 2 D* 


Y - 5 + A) 


'=★1 


d ;2 


TT-" A| 




el2 


+epi 


3D* 


^ + Q*(^) + 

L e*i e*i 


where the subscript “*” denotes evaluation carried out 
at the horizon crossing, and 

Q*(*) = g(^) + (a + 2)QL1ln5 + ^fczH, 

= 2 qL 1 In ^ 

z lo 


(3.38) 


Now we turn to consider the scalar spectral index rig, 
which can be expressed as 

ris — 1 — 2e*i — e *2 ~ 2eli — (3 + 2D*)e*ie*2 ~ D*e*2e*3 


+epi 


3D* 


/C 


*(s) 

-1 


e*i 


AC, 


r(s) 


/cf^' 

e*i 


(3.39) 


where 




--(s) 


The running of the scalar spectral index reads 
as — —2e*ie*2 — e*2e*3 


(3.40) 


+epi 


3D* 


^*(s) 

e*i 


-£ 


*(s) 


.,*(«) 
-I f 


^★1 


(3.41) 


where 


^★(s) 

•'-I 

^★(s) 


= £ 


-i5 


= a£L1ln-+4^\ 


Al(k) c. At{ 1-2(1 +D;)e.^ 


+ + Y— 7 + Aj + 2A5 J e*ie*2 

/ TT^ D*^ \ 

— (4D* + 6 ) (5*He*i + ( Y- 2~ 1 £*2^*3 


+ (2D;^ + 2D; + --5 + At)e*i 


+ 1 “Tip — 2D* + Y ~ ^ ) £*1^*2 


+ epil^) 


(3.43) 


where 

find 


= Qo^- For the tensor spectral index, we 


(3.37) 


Tif ~ —2e*i — 2(1 + D*)e*ie *2 — 2eJ 


+epi 


/3D* 


V 2 


K 


*(t) 


(3.44) 


where ■ For the running of the tensor spectral 

index, we have 


at — — 2 e*ie *2 + epi 


3D* 

2 


,*(t) 
-0 ’ 


(3.45) 


where = 4 *^- 

Finally, with both scalar and tensor spectra given 
above, we can evaluate the tensor-to-scalar ratio at the 
horizon crossing ?;*, and we find that 


16e*i ^ 1 + Dpe *2 + epi 


3D* 

2 


e*i 


(3.46) 


Some remarks about the spectra with the inverse- 
volume corrections now are in order. First, similar to the 
discussions given in [53], as epi oc k~^, both scalar and 
tensor spectra exhibit a deviation from the usual shape 
when k is small enough, i.e., at large scales. Second, as 
a result of the above features in the spectra, the spectral 
indices, and especially the running of the spectral indices 
could be dominated by the quantum gravitational correc¬ 
tions at large scales. Such an interesting feature signals a 
qualitative departure from the inflation given in general 
relativity, and could be crucially important for further 
observational tests of constraints on quantum gravita¬ 
tional corrections. 

Note that in |23l|24| the observables Ug, rit and r were 
calculated up to the first order of the slow-roll param¬ 
eters. Comparing their results with ours, we find that 
they are different. This is mainly due to the following: 
(a) In |23| the horizon crossing was taken as k = H. 
However, due to the quantum gravitational effects, the 
dispersion relation is modified to the forms (Eq. (3.6)) 
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and (Eq. (3.10)), so the horizon crossing should be at 


Wk ='H- Ibl In [23] the mode function was first obtained 
at two limits, k ^ H and k and then matched to¬ 

gether at the horizon crossing where fc ~ "H. This may 
lead to huge errors [31], as neither nor is a 

good approximation of the mode function /i^ at the hori¬ 
zon crossing. This is further supported by considering 
the exact solution for cr = 2 , to be presented in the next 
section. 


IV. INVERSE-VOLUME CORRECTION FOR 
cr = 2: EXACT SOLUTION 


(4.7) 

where 7 = r( 1 -|-zc)/T(-|- ). Comparing the above 

expressions with the initial condition and using the rela¬ 
tion, 

(4.8) 

one gets C 2 = 0 , and 

g 8^02 


Cl = 


y 2 fca'/"' 


(4.9) 


When cr = 2, if we assume that all the slow-roll pa¬ 


rameters are constants, the equations of motion (3.6) and 


(3.10) for both the scalar and tensor perturbations can 


be casted into the form 


Oo 


+ ( ^ + + « 2 fcV ) /ife(^) = 0, (4.1) 


One also can use the Wronskian condition 

k‘k{v)k-k{v) - f^kiv)fJ'kiv) = i 

to determine the coefficient ci, which exactly gives the 
same result. Thus, with the initial condition (4.4) the 
solution reads 


where 


ao = --iy^, ai = 1 - mepiK, 02 = (4-2) 


With Oo, oi, and 02 being constant, Eq. (4.1) can be 
solved exactly, whose solution reads 


k-kiv) = 


e / icLi V . ! _ ,2 2 


yj—2r\k\ 


,1/4 


"''U-4y=-2-v^^V 


( 4 . 10 ) 

Considering the asymptotic form in the limit ry —>■ 0“, 


Mfc(i7) = 


, C 2 


WW - 


-.WM 


4Vdi’ 2 

/ iai V 




^=77 V4v/^’ 2 


1 7^,-iV^k Tj (4.3) 


where WW{hi^h 2 ,z) and WM{hi,h 2 ,z) are the Whit- 
takerW and WhittakerM functions, respectively. To de¬ 
termine the coefficients ci and C 2 , let us consider the 
adiabatic initial condition. 


lim /ife(? 7 ) = , ^ 

-kv^+^ y'2uJk{v) 

1 


ka2^‘^\/—2ri 


,-i j u>k{'n)dri 


^—. (4.4) 


The asymptotic forms of Whittaker functions 
WW{bi,b 2 ,z) and WM( 6 i, 62 , z) in the limit \z\ —>• -l-oo, 
on the other hand, are. 


WW{bi,b 2 ,z) ~ 

WM{bi,b 2 ,z) ~ ^r(l + 262 ) 


r(| +b 2 + bi) 


(4.5) 

(4.6) 


Then, one arrives at 


Cl 


^xk{v) - " i-i^2kW)^ 


r(^) 


p/ 14+1 I iai \ 

^12 4- 4 ^; 


(-zV^fc^r;")^, 


(4.11) 


we find 


k'kiji) ~ 


e 8+02 


r(r 




X [y/^k'^r]' 


2 2\ 


(4.12) 


In the above we have ignored the irrelevant phase factor. 
At this position, the power spectra can be computed in 
the limit 77 —?► 0 “ as 


k^ 


k-k{v) 


z(r]) 


(4.13) 




For scalar perturbations, we find 


Af(fc) = 


Ibir^Mpiei 
2 

‘GR 


~ Aqr ( 1 -I- lymnepi — 


2v{v^ - 1 ) 


-X«^epi 


(4.14) 
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In the above, we have used the asymptotic formula of 


the Gamma function 




ly+l iai ) 


Then consider- 


4^/02 J 

ing the slow-roll expansion of v, Tn(r]), and k for scalar 
perturbation, 


3 e-2 

(4.15) 

ao , o *^ 0^2 

+ 3Q!o r, ; 

ei 2 ei 

(4.16) 

M‘^(I- 2 ei), 

(4.17) 


one finds 


1 + aoH^ 


3 1 

2 e^ 


5 

12 


1 £2 
4 ei 


epi 


(4.18) 


It is easy to show that the coefficient of the leading order 
term epi/ei is exactly consistent with the result obtained 
by the uniform asymptotic approximation, presented in 
the last section. 

Let us turn to consider the tensor perturbations, for 
which we have 




3-2i^ 


_ 1 ) 

X 1 1 + UtTritKepi - - - . Qg/^epi 


(4.19) 


Considering the slow-roll expansions of vt iv) j (v) > 

K, 


3 

- ao + aoci. 


(4.20) 

(4.21) 


we find 

A2(fc) ~ [1 - aoHhpi - SaoHhiepi] . (4.22) 

One can check that the coefficient of the leading order 
term from the uniform asymptotic approximation ob¬ 
tained in the last section is — ^ — 1 . 011 , which is very 

close to the exact value — 1 , obtained from the above ex¬ 
act solution. Therefore, the results presented in the last 
section for cr = 2 are the same as these exact results ob¬ 
tained in this section within the errors allowed by our 
approximations. 


V. DETECTABILITY OF QUANTUM 
GRAVITATIONAL EFFECTS 

With the slow-roll conditions, the holonomy correc¬ 
tions are normally much weaker than the inverse-volume 
ones, and their effects in the early universe are not ex¬ 
pected to be observed in the near future m- However, 


this may not be the case for the inverse-volume correc¬ 
tions [32]. Therefore, in this section we shall consider 
only the latter. 

In general, with the spectral index and running given 
in Eqs.(3.39) and (3.41), the scalar spectrum can be ex¬ 
panded about a pivot scale ko as 


InA^(fc) ~ lnA^(fco) + [ns(fco)-l]ln — 

fco 


s{ko) 


fco 


oo (n) 

E O^s 

Ti 

n=3 


(fco) 


In" ( f 


(5.1) 


where up to the second-order approximations in terms of 
the slow-roll parameters and the leading order contribu¬ 
tion from epi we have 


n-s — 1 — —2e*i — e*2 ~ ~ (3 + 2fcln) e*ie*2 

—u„e*2e*3 + epi I j Gi i 


C4s — —2e*ie*2 — £*26*3 ~ o'epi 


SET* 


„*(s) -1 
^-1 £*1 ! 


(5.2) 


and 


a(")(fc) (-l)"-V"-iep, (5.3) 

Note that when a — 3, one has to replace by 

Similar to [23], it is easy to find 

cts \ko) ^ 


m—S 


n: 


= —epi 


fco 

3H* 

2 


^*(s) 1 

'^-i £*i 


/, fc , 2 - 

X In --— In — 

\ Kq 2 Kq 


k e '^0 — 1 


)(5.4) 


In order to carry out the CMB likelihood analysis, it is 
convenient to introduce the following potential slow-roll 
parameters. 


Ml; U| 


MpiV^v .2 _ 

) Kv = 


V 


V2 


(5.5) 


with which the scalar spectrum can be cast in the form 
A^(fc) ~ A^(fco)exp 


/ X , fc , a* fco , 2 fc 

{n,(ko) - l)ln — + — —In — 
fco 2 fco 


2 <T -i 


ev 

, fc cr 2 k ^ - 1 

X In --- In — H- 
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(5.6) 

while the spectral index and its running can be written 
in the forms, 


TT-s 1. In the forthcoming experiments, specially the 
Stage IV ones, the errors of the measurements on both 
Ug and r are cr(ns), (y{r) < 10“^ |39j . which implies the 
error of the measurement of r„(ns,r) is 


~ 1 — 6ev + 2r]v — 


10 

T 


24D( 




- 2)evVv + ^+{--2D:]C 


ev 


JC 


2 ri 

3 

*(s) , cr2((T-3)ao 
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mcT-a-3) 

(5.7) 


and 


cts — ~‘2 ^y + 16?7yey — 24ey 
aepiHl 


ev 


(5.8) 

2<T -1 0 ( y J 


Again remember that when cr = 3 one has to replace 
/C_Ve*i by Ap' h 

Let us consider the power-law potential 

V(^) = Voyp^, (5.9) 

where Vq and n are constants. In this case, it follows that 


ey = 


_ 


rjv = Mpi- 


— -^Pl 


4 n^(n — l)(n — 2) 




(5.10) 


from which one can reduce the potential slow-roll param¬ 
eters to one (i.e., ey). 


2(n-l) 

Vv = -ey, ?y = 


2 _ 4(n- l)(n- 2) 2 

„2 ^V- 


(5.11) 


It is also convenient to parameterize the inverse-volume 
corrections as 


6{k) = aoepiH'^. 


(5.12) 


Thus, in the scalar power spectrum, with the power-law 
potential there are only two independent free parameters, 
ev{ko) and (5(fco). 

Now let us turn to consider the observational effects 
of the inverse-volume corrections. When inverse-volume 
contributions vanish, we have rig = ng(ev) and r = r{ev)- 
Thus it is easy to show that, up to the second-order of 
ey, the relation [57], 


r„(ns,r) = (ug - 1) -I- 


(2 -I- n)r 
8n 


(3n" + 18n-4)(n,-l)\ 

^ 6(n + 2)2 

holds precisely. The results from Planck 2015 are Ug = 
0.968±0.006 and ro,oo 2 < 0.11(95% CL) [55] . which yields 


a(r„) < 10 


-3 


(5.14) 


Therefore, if any corrections to Ug and r lead to 
rn(?^s, r) > 10“^, they should be within the range of de¬ 
tection of the current and forthcoming observations [39] . 

In particular, when the inverse-volume corrections are 
taken into account (Jpi ^ 0), w e have Ug = ns(ey,epi) 
and r = r(ey, epi), and Eq.(|5.13 1 is modified to. 


Tn{ng,r) = F{a) 


m 


ey 


(5.15) 


where 5{k) = aoCpiTL®’ and 


J-(a) = 


(T^(cr — 3)ao 
18 


(3D>- cr - 3). 


(5.16) 


Clearly, the right-hand side of the above equation repre¬ 
sents the quantum gravitational effects from the inverse- 
volume corrections. If it is equal or greater than 0(10“^), 
these effects shall be within the detection of the current 
or forthcoming experiments. It is interesting to note that 
the quantum gravitational effects are enhanced by a fac¬ 
tor of which is absent in [55]. 

In the following, we run the Cosmological Monte Carlo 
(CosmoMC) code [40] with the Planck [41], BAO [42] . 
and Supernov a Le gacy Survey [35] data for the power- 


law potential (5.9| for n = 1, ^ respectively. It 


is worthwhile to mention that all these potentials can 
be naturally realized in the axion monodromy inflation 
motivated by string/M theory [H]. In [55], by using 
CosmoMC code [30] , we already extracted constraints on 
loop quantum correction parameter dpi /ey and slow-roll 
parameter ey for cr = 1 and a = 2 when n = 1. It was 
noted that the constraint for cr = 2 is much more tighter 
than the case of cr = I, and makes the quantum gravi¬ 
tational effects undetectable for models with cr > 2. It 
is interesting to note that models with small cr are also 
favorable theoretically (55] [53]. Therefore, in the follow¬ 
ing we shall focus on the observational constraints for the 
case cr = I. 

We assume the flat cold dark matter model with effec¬ 
tive number of neutrinos A^e// = 3.046 and fix the total 
neutrino mass = 0.06eld. We vary the seven param¬ 
eters: (i) baryon density parameter, (ii) dark mat¬ 
ter density parameter, (hi) the ratio of the sound 

horiozn to the angular diameter, 0, (iv) the reionization 
optical depth r, (v) d(fco)/ey, (vi) ey, (vii) A‘^{ko). We 
take the pivot wave number ko = 0.05 Mpc“^ used in 
Planck to constrain 6(kg) and ey. In Fig{^ the con¬ 
straints on d/ey and ey are given, respectively, for n = 1, 
n = |, n = |, and n = ^. In particular, we find that at 
68% C.L., 


Hko) < 7.9 X 10 


-5 


with n = 1, 


(5.17) 
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FIG. 1: Two-dimensional marginalized distribution for the parameters 5/ey and ey at the pivot ko = 0.05Mpc“^ for the 
power-law potential with n=l, n=|, n=|, and n = respectively. The internal and external lines correspond to the 
confidence levels of 68% and 95%, respectively. 


S(k(}) < 5.7 X 10 with n = -, (5.18) 

O 

S(k(}) < 5.4 X 10“®, with n = -, (5.19) 

5 

S(ko) ^ 3.4 X 10“®, with n = ^, (5.20) 

o 

which are much tighter than those given in |24j . In the 
above, to get the the bounds on S{ko), we have used 
the best-fit values of ey, respectively for different n. It 
is also easy to see that the upper bound for both (5(fco) 
and ey decrease slightly as n decreases. However, it is 
remarkable to note that the up bound on S{ko)/ey is 
rather robust for different values of n, which is roughly 

~ 0(1) X 10-3 (68% C.L.). (5.21) 

ey 


With such bounds, the gravit ation al quantum effects, de¬ 
noted by J^((t)^^^ in Eq. (5.151, could be within the 


range of the detection of the current and forthcoming 
cosmological experiments |39j . In addition, as we already 
pointed out in [35], the up bounds for increase dra¬ 
matically as cr decreases. Thus it is very promising to 
expect the detectability of gravitational quantum effects 
for cr < 1 in the forthcoming cosmological experiments. 


VI. CONCLUSIONS 

The uniform asymptotic approximation method pro¬ 
vides a powerful, systematically improvable, and error- 
controlled approach to construct accurate analytical so¬ 
lutions of linear perturbations. Its effectiveness has been 
verified by applying it to the inflation models with non¬ 
linear dispersion relations [29] and A:-inflation [30]. In 
this paper, we apply the high-order uniform asymptotic 
approximation to derive the inflationary observables for 
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scalar and tensor perturbations in LQC with holonomy 
and inverse-volume quantum corrections. We obtain 
explicitly the analytical expressions of power spectra, 
spectral indices, and running of spectral indices up to 
the third-order approximation in terms of the parame¬ 
ters introduced in the uniform asymptotic approxima¬ 
tion method. To this order, the upper error bounds are 
< 0.15%, accurate enough for the current and forthcom¬ 
ing experiments |3Hj . These expressions are all descibed 
in terms of the slow-roll parameters (up to the second- 
order) and the parameters which represent the holonomy 
and inverse-volume quantum gravitational corrections. 

For later applications of our results, we also rewrite 
all the inflationary observables including power spectra, 
spectral indices, and running of spectral indices for both 
scalar and tensor perturbations in terms of quantities 
evaluated at the time when the inflationary scalar (or 
tensor mode) crosses the Hubble horizon. With the re¬ 
sulting expressions, the tensor-to-scalar ratio is also ob¬ 
tained, and it is shown that the holonomy corrections 
do not contribute to the tensor-to-scalar ratio up to the 
second-order approximations. More interestingly, with 
the inverse-volume corrections, we find that both scalar 
and tensor spectra exhibit a deviation from the usual 
shape at large scales, which could be potentially impor¬ 
tant for the observational tests. As the uniform asymp¬ 
totic approximate solution at the third-order has error 
bounds < 0.15%, the inflationary observables obtained 
in the present paper represent the most accurate results 
obtained so far in the literature. 

Utilizing the most accurate CMB, BAO and SN data 
currently available publicly [iTHiH] , we also carry out the 
CMB likelihood analysis, and find the tightest constraints 
on ((5(/co)) ey)) obtained so far in the literature. Even 
with such tight constraints, the quantum gravitational 
effects due to the inverse-volume corrections of LQC can 
be within the range of the detection of the current and 
forthcoming cosmological experiments |39) , provided that 
cr < 1. 
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Appendix A: The uniform asymptotic approximation 

In this section, we present a brief introduction of the 
uniform asymptotic approximation method and its appli¬ 
cations to the inflationary cosmology for the cases where 


the dispersion relation has only a single turning point. 
Eor details, we refer readers to our original papers [23- 

ISOl- 

A. The approximate solution of the mode function 

In the uniform asymptotic approximation method 
[2312811311 Eg, one usually works with the second-order 
differential equation, 

+ 9(y)] Mv)- (A.i) 

In the above the parameter A is used to trace the order of 
the uniform approximations. Usually A is supposed to be 
large, and it also can be absorbed into g{y). Thus when 
we turn to determine the final results, we can set A = I 
for the sake of simplicity. Eor convenience, we also use 
the notation g{y) = \^g{y). Specific to the cosmologi¬ 
cal applications, Pk{y) represents the inflationary mode 
function for cosmological scalar or tensor perturbations, 
and one can identify 

A'5(y) + q{y) = -^ , (A.2) 

where is the associated dispersion relation for the 

inflationary mode function pLk{y) and z{ri) depends on 
the cosmological background evolution. In most of the 
cases, g{y) and q{y) have two poles (singularities): one is 
at 2 / = O'*" and the other is at y = -|-oo. As we discussed 
in [22] (see also [33 Eg), if these two poles are both 
second-order or higher, one needs to choose 

9(2/) = -^’ (A-3) 

for ensuring the convergence of the error control func¬ 
tions. In this paper we shall restrict our investigations 
to this case. In addition, the function g{y) can vanish 
at various points, which are called turning points or ze¬ 
ros, and the approximate solution of the mode function 
/ifc(y) depends on the behavior of g{y) and q{y) near these 
turning points. 

To proceed further, let us first introduce the Liouville 
transformation with two new variables U{^) and ^ via 
the relations, 

u{f) = (A.4) 

where y = ^' = d^/dy, and 

fiO = fV\My)\dy, /W(0 = ^. (A.5) 

Note that x must be regular and not vanish in the in¬ 
tervals of interest. Consequently, /(^) must be chosen so 
that has zeros and singularities of the same type 
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as that of g{y). As shown below, such a requirement 
plays an essential role in determining the approximate 
solutions. In terms of U and Eq. (A.l) takes the form 


cPU 

~d^ 


u, 


where 


'0(0 = 


X dy"^ 


(A.6) 


(A.7) 


and the signs “±” correspond to g{y) > 0 and g{y) < 0, 
respectively. Considering = 0 as the first-order ap¬ 
proximation, one can choose so that the first-order 

approximation can be as close to the exact solution as 
possible with the guidelines of the error functions con¬ 
structed below, and then solve it in terms of known func¬ 
tions. Clearly, such a choice sensitively depends on the 
behavior of the functions g[y) and q{y) near the poles 
and turning points. 

In this paper, we consider only the case in which g{y) 
has only one single turning point j/o (for g{y) having 
several different turning points or one multiple-turning 
point, see [35]), ke., g{yo) = 0. In this case we can choose 




(A.8) 


where ^ = ^{y) is a monotone decreasing function, and 
± correspond to g{y) > 0 and g{y) < 0, respect ively . 
Following Olver the general solution of Eq. (A.6) 
can be written as 


= CIO 


Ai(A='y)y: 


A.d) 


s=0 

Ai'(A2/3e) 


A2« 


E 


+/3o 


A4/3 

n 

s^O 

, Bi'(A2/3^)g 

s=0 


s=0 

M) 

A2s 


BsiO 

X2s 


(2n+l) 


A4/3 


Ml 

A2s 


+ e 


(2n+l) 


where Ai(a;) and Bi(a;) represent the Airy functions, 
^( 2 ra+i) errors of the approximate solu¬ 

tion, and 

Ao(0) = 1, 

BsiO = 




As-hi(0 = -\Bs{ 0 + ^ J ^{v)Bs{v)dv, 


(A.IO) 


where ± correspond to ^ > 0 and ^ < 0, respectively. 
The error bounds of and can be expressed 


as 


(2n-|-l) 


de 


(2n+l) 


/5C 


M(A2/3^)’ A2/3Ar(A2/3^) 


< 2E-i(AMexp 






A'2n+1) 


de 


X2n+1 

(2n+l) 


/d^ 


M(A2/3^)’ A2/3iV(A2/3^) 


< 2E{\^/^^) exp 


2Koy^A\e'^\Bo) 




X2n+1 


(A.Il) 


where the definitions of M{x), N{x), kq, and 'f'a,b{x) can 
be found in [28l. 


B. Power spectra and spectral indices up to the 
third-order 


With the approximate solution given in the last sec¬ 
tion, now let us begin to calculate the inflationary power 
spectra and spectral indices from the approximate solu¬ 
tion. We assume that the universe was initially at the 
adiabatic vacuum, 

lim gk{y) = lim d^ (a.I2) 

y^+ca ^2uJkig) 

Then, we need to match this initial state with the approx¬ 


imate solution (A.9). However, the approximate solution 


(A.9| involves many high-order terms, which are compli¬ 


cated and not easy to handle. In order to simplify them, 
we first study their behavior in the limit y —>■ -|-oo. Let us 
start with the Bo{^) term in Eq. ( A.lOj), which satisfies 


BoiO = - 


1 


2Mio M 2M 


(A.I3) 


where dv'il^{v)/\v\^y is the associ ated error 


(A.9) control function of the approximate solution (A.9), and 


in the above we have used Ao(^) = 1. The error control 
function A^(^) is well behaved around the turning point 
yo and converges when y —>■ -|-oo. As a result, we have 


hin BoiO = 
y-y-l-oo 2y—4 

Then, let us turn to Ai, which is 

rC. 


(A.I4) 


Ai{0 = -\Bo{ 0 + \ '>p{v)Bo{v)dv. (A.I5) 
In the limit y —>■ -l-oo, Bq(^) vanishes, and we find 


I r 

lim Ai(^) = -X / 
/^+oo 2 Jq 


li 

r j 

1 

> 

o 

1 

> 

o 


dv 
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J^(—oo) 


(A.16) 


Note that in the above we have used the formula 

/•?2 


n! t f {in) f{in-l)' 


t 

'io 


f{ii)diidi2---di„ 


f{v)dv 


(A.17) 


Thus, up to the third-order, we have 


oU) + ^2 2A2 


J^{—oo) 


+ 0 - 


Boii) 

X 


1 J^{—oo) 

2A 


O 


A3 


Using the asymptotic form of Airy functions in the 
limit ^—>■—00, and comparing the solution fj,k{y) with 
the initial state, we obtain 


ao = 

Po = i 


Ai/3 

2k (Ao + Ai/A2)-zy=eSo/A’ 

TT A3/3 

2k (Ao + Ai/A2)-zv/^Bo/A’ 

where we have 

(Ao + Ai/A^) - iy^Bo/X = (1 + 0(1/A3))e*^ 


(A.19) 


(A.20) 


Here 9 is an irrelevant phase factor, and without loss of 
generality, we can set 0 = 0. Thus, we finally get 


ao 

Ai73 


Po 

2k' Ai/3 


TT 


(A.21) 


After determining the coefficients ao and /?o, we can 
calculate the power spectra of the perturbations. As j/ —)■ 
0, only the growing mode is relevant. Thus we have 


Mfc(j/) - 




9{y) 


1/4 




s=0 

+00 


A2/3Bi'(A2/3^) ^ Bsii) 

2^ X2s 
s=0 


A2 


. (A.22) 


In order to calculate the power spectra to higher order, 
let us first consider the Bo{i) term, which satisfies 


1 


(A.23) 


2 (1/2 2(1/2 

In the above we had used the relation ip^'^di = —y/gdy. 
Knowing the Bq term, we can get the Ai term, which is 

A 1 'P{v) r P{u) J 

= ii 




(A.24) 


Thus up to the third order and considering the asymp¬ 
totic forms of the Airy functions in the limit i —)■ -l-oo, 
we find 


Poei^^ 

= Xl/Ogl/V/2 


J^{+oo) J^{+oo)'^ 


2X 8A2 

C>(1/A3)|. (A.25) 


Then, the power spectra can be calculated, and is given 

by 


. (A.18) ^ 


k^ 

27r2 


yk{y) 


l/-i-0 + 


1 +A<±^ + ^d±2d+o(i/A3) 


2A2 


(A.26) 


From the power spectra presented above, one can get the 
general expression of the spectral indices, which now is 
given by 


zz — 1 = 


din A2(fc) 


dlnfc 


~ 3-b2 


din A: 


^-Tx- 1 d,^(-|-oo) 


+<31^ 


(A.27) 


where the last term in the above expression represents 
the second- and third-order approximations. It should 
be noted that the above results represent the most gen¬ 
eral expressions of the power spectra and spectral indices 
of perturbations for the case that has only one-turning 
point. 


Appendix B: Integral of \fg{y) and the error control 
fnnction M’{-\-oo) with inverse-volume corrections 

In general, the integral of ^/g can be divided into two 
parts, 

fVo 

/ V9{y)dy- Ii + I 2 , (B.l) 

where after some tedious calculations we find 


lim/i = yo 

jz-i-o 


-1-ln^-f (-^-fln2-ln:^ )Ao 
2yo V 2 yo 


+ l(-2 + 7r)Ai + (l-0 A 2 
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+ 

+ 

+ 


2 7r\ 2 SttA ^ 

8 37r\ , / 8 57r\ , 

■l5 + l6)^^+(l5"^)^® 


16 57 r\ 

-^— M? 

35 32 y 


(B. 2 ) 


Br = -^{e-f)€piyQ ^ (ko™! + kiTOq) . 
Bs = -^/epiyo”^ (^owi + KiTOo) • 


(B. 5 ) 


The error control function ,^(+00) can also be ob¬ 
tained by employing the expansions given in Eq. ( 3 . 271 , 
and we Hnd 


lim I 2 = — — 


y-J-O 


24 


+2/0 


In 2 1 2 2 / , - 

^5 -^ In — 2/1 

2 2 yo/ 


TT” , „ In 2 1,22/ 


— — — In 2 
24 


2 

TT TT In 2 \ „ 

-+ - ^3 

2 2 / ^ 


-i?i — 7^2 In 2 

2 ^ ^ 

, , , f biT 37 r In 2 \ 

-|-(1 - ln 4 ) i ?4 + ”1 - ^ J ^5 

14 81 n 2 \ /I 57 rln 2 477 r\ 


148 16 In 2 


B. 


75 5 

where +0, • • ■ , Bq, • • • , Bg, and Cq are given by 


(B. 3 ) 


1 


n ) — 2 - 1-(7 


-amo + xvl) epiKo, 


^0 = 2% 

^1 = (-^TOo + xvl) epi^o, 


— 2 - 1-(7 


(-CTOo + axvl) epiKo, 


^2 = 2% 

-43 = 2^oepiyo”^ (^X2/o “ ^"^ 0 ), 

-44 = ^Koepiyo”^ (cxyo “ e^o), 

-45 = ^Koepi2/o”^ - /^o), 
. 1 

-46 = ij^eKQX^pWQ , 

- 4 ? = ^/Koxepiyo. 


(B. 4 ) 


and 


50 = --epiyp" («odT.i + Kirho), 

51 = ^(a - 5)epiyo“^ (kot+i-I-/ cimo), 

52 = ^epiyo“^(a-f 6 - c) (ko+ii- f Ki? 7 io), 

5 3 = -^epiyo“^(Ko77ii-f Ki?Tio)(a-6-c-f d), 

5 4 = -^epiyo“^(KoTOi+KiWo)(^-c-d + e), 

5 5 = -^epiyg”^ (koWi + Ki+io) (c-d-e-I-/), 
Be = -iepiyo“^(d-e -/) (Ko+ii + Ki?Tio), 


=^(-|-oo) 

A 


1 ( 23 -f 121 n 2 )yi 

6^0 


72%" 
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Tra 

("2 

TT \ / 

—— -|- 

- - 

757 ^ + 

16 

V 3 

16 / \ 
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3 ~ 
32 
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IStt 
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d H" 


/17571 8 


- X e 


1757 r\ 1 

128 ) ^ 12 


Xepi2/o 


(a 2d 2e 8/ ctt dir 

+-^- - -^- 

ll 2 3 3 3 16 16 


15e7r 15/7r \ _ _ 

1 _ I - rt I ^ 
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32 

23 

lU ^ 

1 

12 
7tt 
48 


32 




woepiKodo 


TT 


In 2 

U 

7 tt ^ 7rln2 
48 16 

1 2 In 2 

6 lY ' 


1 

12 


-V 

48 J 


1 4l7r 

6 ~ lY 

4 4l7r 
9 lY 


TT In 2 
16 

21 n 2 157 rln 2 


3 

81 n 2 


32 

157 r In 2 


3 32 

X (koTOi -f KiT+o) epiyY^- 


(B. 6 ) 


Once we get the integral of \ /g{y) in Eq. (B. 2 ) and the 
error control function in Eq. (B. 6 ), from Eq. (A. 261 we 


can easily calculate the power spectra. 

Now we turn to consider the corresponding spectral 
indices. In order to do this, we hrst specify the k- 
dependence of yoivo), ViiVo), epi(do), and fh{yo) through 
do = do(fc)- From g{yo) = 0 we hnd 

% - % + ^epi(wof'o~^ “ (B.7) 

and noticing —krjo = yo{r]o) and epi ^ we obtain 


dln(-77o)_f 1 + ^ 


dlnfc 


2/0 


1+2^pi(’^oY ^-XPo)«^o 


(B. 8 ) 
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Then, using the above relation, the spectral indices are 
given by 

n — 1 ~ 3 — 2i>o + ( — In 4 I Di 

VK J 

+epi (miKo + ffioKi) + £ 21^0 

+epiX^o (^3^0 ^ 

+epiXl^i (^5^0 ^ 

+epimoKo (Syi/g ^ + Sgi^o > (B-9) 

where Si, i = 1 , • • • 8 depend on the value of tr and are 
given in the Table I. Then the corresponding spectral 
index reads 

a ~ 2 j>i - crepimoKo (Pq~^Et + “^Eg) 

-epi(miKo + TOoKi)[t'o~^ (o’Si + E7) 

+ ^0 ^ (o '^2 + Sg) 

-CTxepiKo + *^ 0 ^^ 514 ) . (B.IO) 


+12e2eiSH — 2,e2^i5H — ‘2£2£3 ^iSh + 7e2ei 
11 


+4e2e3e4 + —e2€i 


17 

+ -^e2e3e4ei 

0 


+ 


e2e3e4e5 


17 2 .2 

"^e2e3ei + 6e2e3ei 


e2e3 


£2 £3 £4 


£2£3£4 


(C.3) 


and 


s 3 q 2 e3e2 2 e 3 £ 4£2 

^3 — ~£i £2 ~ 3£i£3£2- 2 -£i*^3£2- 2 — 


:£3£4£2 


1 


£i£3£4£2 


:£3£4£5£2- 


(C.4) 


For the tensor perturbations, using Eq. (2.12), we find 

+ 4 +- eSj^ei + el 
2^2^^ 34 2 4 2 4 

— ^07464 — - 0 _f/£l £2 + + 2 ^ 1^2 + 2£i£2£3, 

(C.5) 


Appendix C: Slow-roll expansions of v(r\), Ca{ri), m{rj), 
K{ri), and their derivatives 

A. Expansions with the holonomy corrections 


Let us first consider the scalar perturbation with the 
holonomy corrections. Using the expression of z"/z in 
Eq. (2.111, it is easy to find that 1 /( 77 ) foi' scalar pertur¬ 


bations reads 

S/ \ ^1 I ^2 2 I ll£2£l £2£3 

i/ ( 77 ) ~ 2 + 


2el5H 


17 


— QeiS^j — 

14 


+ -^£ 2 £i + -j^£2£3£i - 


6 

2£2£i^/4 

3 

£^£3 


6 


— 2eiSH 

77 

"18 


£ 2 £? 


18 


(C.l) 


Consideration of the derivatives of 12 ^ [t]) with respect to 
In (— 77 ) yields 


and 


—£i£2 — 'iSnel + 2(5preie2 ~ ^4^2 — ^^ 1^2 

~2£i£2£3 + 286ljel H —~ 

-25Hele2 + be ?£2 - ^(5ff£i£2 + y£?£2 
432 ^ 46 2 ,2 

4 94 

-\--^^1^2^3 + 2ei^2e3e4, (C. 6 ) 

£i£2 + £i£2£3 ~ SSnel + 12SHele2 — 25Heie2 
+7£4£2 + 2 ^i ^2 ~ 2(5preie2£3 -I- 4e4e2£3 

+4£i£2£3 + 2£i£ 2£3 + g£l£2£3£4, (C.7) 
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(C.2) 


— —£i£2 ~ 3616362 — 616362 — 61636462. 


(C. 8 ) 


Now we turn to consider 03 ( 77 ). Expanding it in terms 
of 5h, we observe that 


63 ( 77 ) = a/i - 25H ~1-5h - Y - 

Similar to 12 ^ (rj), the derivatives of £ 3 ( 77 ) with respect to 
In (— 77 ) are given by 

61 ( 77 ) ~ -261(5/7 - 26i5/7 - 46i5|i - 96i(5|i - 46i(5|i 

-2ele2SH - 2elSH, (C.IO) 
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TABLE I: Values of Coefficients Ei(i = 1,..., 8) for different values of a. 
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+ 12e2£iSH — — 2e2€3iiS[{ ■ 

(C.12) 


Note that at the turning point we write m{rio) = fho- 
Now consideration of the derivatives of rn{r]) with respect 
to ln(— 77 ) yields 


Note that for both scalar and tensor perturbations, the ^ (T^ao(3 — cr)e 2 

effective sound speed Cs{r]) takes the same form. Thus in ' 

the following, we are not going to distinguish them. 

Similarly, for tensor perturbations, we have 


(C.14) 


B. Expansions with inverse-volume corrections 


Now let us turn to consider the slow-roll expansions 
of ly, m(ri) and K{ii). For v(r]) and its derivatives, it 
is worth to note that when one considers the inverse- 
volume corrections, i'{ri) and its derivatives with respect 
to ln(— 77 ), i.e., vq, V 2 , and V 3 all take the same form 
as those given in general relativity. Thus, they can be 
directly obtained from Eqs. ( C.l l-( C.4 ) for the scalar 
perturbations and from Eqs. (C.5)-(C.8| for the tensor 


perturbations by taking the holonomy parameter 6 h = 
0. For the scalar perturbations, on the other hand, the 
function m(rf) reads 


= 


(3 — a) Oq 
4ei 


3(7 



2 

4 

12 j 


, 5o-^ 


ao 


m iv) - -(C.15) 

while up to the second-order in the slow-roll parameters 
we have ml — 0 and ~ 0. Note that we also write 
m*(rio) as ttiq. 

To get the slow-roll expansions of the inflationary ob¬ 
servables, we also need the slow-roll expansions of K{r]) 
and its derivatives, which are given by 

K{r]) = (C.16) 


Ki = 


(C.17) 


Note that at the turning point we write ^( 770 ) = kq- 
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